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Abstract 

We consider the model of random sequential adsorption, with depositing objects, as well as 
those already at the surface, decreasing in size according to a specified time dependence, from a 
larger initial value to a finite value in the large time limit. Numerical Monte Carlo simulations of 
two-dimensional deposition of disks and one-dimensional deposition of segments are reported for 
the density-density correlation function and gap-size distribution function, respectively. Analyti- 
cal considerations supplement numerical results in the one-dimensional case. We investigate the 
correlation hole — the depletion of correlation functions near contact and, for the present model, 
their vanishing at contact — that opens up at finite times, as well as its closing and reemergence 
of the logarithmic divergence of correlation properties at contact in the large time limit. 

PACS numbers: 68.43.Mn, 02.50.~r, 05.10.Ln, 05.70.Ln 
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I. INTRODUCTION 



Random sequential adsorption (RSA) model has attracted a lot of attention and has a long 
history l|, Q, 3, 3- It finds applications, e.g., fl, jj, fl, IJ, fl, J 3, 9] in many fields, 
ranging from surface science to polymers, biology, device physics, and physical chemistry. 
Traditionally, the RSA model assumed that particles are transported to a substrate which 
is a continuous surface or a lattice, the latter convenient for numerical simulations. Upon 
arrival at the surface the particles are irreversibly deposited, but only provided they do not 
overlap previously deposited objects. Otherwise, the deposition attempt is rejected and the 
arriving particle is assumed transported away from the substrate. The original RSA model 
studies were largely motivated by surface deposition of micron-size objects, such as colloid 
particles. Particles of this size are typically not equilibrated on the surface, and are larger 
than most surface features and the range of most of the particle-particle and particle-surface 
interactions. Various generalizations of the basic RSA model have been considered in the 
iterature, e.g . 
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Q, 0, Q, H, Q, Q, Q, Oh H, Q, Q, Q, Q, Q» Q» 

29j, including "soft" rather than hard-core particle-particle interactions, as 
well as relaxation by motion of particles on the surface. 

Recently, experimental surface-deposition work has expanded to nano-size particles and 
sub-micron-feature patterned (ultimately, nano-patterned) surfaces {^(j], [31], 32], 33]. The 
added control of the particle and surface "preparation" as part of the deposition process could 
allow new functionalities in applications, and therefore it has prompted new research efforts. 
Specifically, deposition on surfaces prepared with patterns other than regular lattices was 



studied 24|, 25| motivated by new experimental capabilities in surface patterning. Another 
development involved a study 8], [lsj] , [l?]] of one-dimensional deposition of segments that, 
after attachment to the substrate, can shrink or expand, motivated by potential applications, 
e.g., in device physics [n]], 

The motivation for our present work has been a newly emerging experimental capability 



34j, 35[ of depositing polymer "blobs" (or polymer-coated particles) the size of which can 



be modified by changing the solution chemistry: Approximately spherical particles can be 
deposited in a process whereby their effective size (including the interaction radius) is varied 
on a time scale comparable to that of the deposit formation. The size of both the particles 
in solution and those already deposited will thus vary with time, in a controllable fashion. 
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For particle deposits formed on pre-patterned surfaces, an interesting property is sug - 
gested by experiments, 3(|, {38], (^J, 40], 41], 42], and theoretically verified, 24], 



251 ]: They acquire semi-ordering properties "imprinted" by the substrate, as quantified by 



the development of peaks in two-particle correlation. However, as in the original irreversible 
RSA model, there remains a significant peak at particle-particle contact, which for depo- 
sition continuing indefinitely, at infinite times becomes a weak singularity 43j, 44J. This 
tendency of particles to form clumps, is undesirable mostly because many nanotechnology 
applications rely on nanoparticles utilized in isolation, or simply being kept away from each 
other to avoid merging. In this work, we establish that deposition of particles of varying, 
specifically, decreasing size, can yield deposits without clumping, by opening a "correlation 
hole," i.e., a property of depletion of two-particle correlations near contact, as defined in 



other fields, e.g., 



45|, |46j, j47j, [48j, [49(, [50|. In our case, the correlation functions studied 



actually vanish at contact for finite deposition times. 

The outline of this article is as follows. In Section[TTl we consider the two-dimensional (2D) 
deposition of shrinking disks in a plane. The model is defined and then a two-point density- 
density correlation function is studied by numerical Monte Carlo simulations. In Section [TTTl 
we present both numerical and analytical analysis of the one-dimensional (ID) deposition of 
shrinking segments on a line. Our ID study focuses on the gap-density distribution function. 
A brief summary is offered at the end of Section IIHI 



II. DEPOSITION OF DISKS ON A TWO-DIMENSIONAL SUBSTRATE 



A. Definition of the model 



In this section, we introduce the model for the most relevant geometry for possible appli- 
cations: The 2D RSA model of deposition of disks of diameters D(t) on an initially empty 
planar substrate. A ID model of deposition of segments on a line provides additional insight 
into the problem and will be considered in the next section. 

As usual in RSA, we assume that disks are transported to the surface with the resulting 
deposition attempt flux R, per unit area and unit time, t. A disk adsorbs only if it does 
not overlap a previously deposited one. The diameter of all the disks, those already on the 
surface and those arriving, is a decreasing function of time, D(t) < 0, varying between two 
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FIG. 1: The left panel illustrates deposition attempts in 2D and ID that are rejected due to particle 
overlap. The right panel shows a situation at a later time at which the same deposition attempts 
would succeed due to particle shrinkage. 

nonzero values -D(O) > D(oo) > 0. For the 2D model, we carried out numerical Monte 
Carlo simulations to estimate the particle-particle pair correlation function that describes 
the relative positioning of the disks with respect to each other. It is defined as the ratio of 
the number of particle centers N at distances from r to r + dr from a given particle's center, 
normalized per unit area and per the deposit density, p, 

P( rf )- N(r,dr,t) 

In studies of RSA, p{t) is usually a quantity of interest. It grows linearly for short times 
and reaches a jamming-limit value, which is less than close packing, for large times. In the 



latter regime the 



approac 



i to the jamming coverage on a continuum substrate is described 



by a power law [43[, [44j], [5l|- However, in our present study we found by numerical 
simulations no new interesting features for pit), for the considered time- dependences of the 
disk diameters (specified below). Therefore, we focus our presentation on the correlation 
function, and specifically its properties near particle contact, at r ~ D(t). 

The model is illustrated in Figured) We note that if the disks shrink too fast as compared 
to the time scale of the buildup of the deposited layer, which is of order ~ 1/RD 2 , then 
the problem will be reduced to that of simply depositing smaller disks. On the other hand, 
if the disks shrink too slow as compared to the time scale ~ 1/RD 2 , then the depletion 
of the correlations at contact will be trivially attributable to disk shrinkage alone. We are 



interested in the interplay of two effects: The decreasing disk size enlarges voids between 
already deposited disks. At the same time this process increases the rate of successful disk 
deposition events, which reduce voids between disks. Given that numerical simulations for 
this problem are quite demanding, we report results for the following time-dependence for 
the disk radius, 

D(t) = L>(oo){l + exp[- J R J D 2 (oc)t]} . (2) 
We also studied numerically the cases 

D(t) = D{ cc){ 1 + WT ±_ }, (3) 

CW = C M{i + i + ^ 2(oo) J , (4) 

and found qualitatively similar results. 
B. Numerical results 

Let us now briefly outline the numerical procedure used in our Monte Carlo simulations. 
In order to calculate the correlation function P 2 (r,t), we have to generate a distribution 
of deposited objects. We used periodic boundary conditions and system sizes L x L = 
500-D(oo) x 500-D(oo). The ^-coordinates of the center of the next disc which makes 
a deposition attempt were randomly generated. If the disk does not overlap any of the 
previously deposited ones, then it is deposited. The total number of deposition attempts is 
thus equal to RL 2 At for a physical time interval At. All our data presented below, were 
averaged over 100 independent Monte Carlo runs. 



We point out that there exist algorithms for studying RSA at large times [52J, specifically 
designed to account for the fact that isolated residual voids (defined as landing areas for 
particle centers for allowed depositions) then evolve independently. However, in our case 
we were primarily interested in intermediate times. Furthermore, as disks already adsorbed 
on the surface shrink, all the voids increase and some might merge to form larger voids. 
Therefore, we used the straightforward algorithm described in the preceding paragraph. 

Our main result is illustrated in Figure El The correlation function introduced earlier, is 
plotted for several times, which are multiples of a characteristic process time-scale 

r = l/RD 2 (oo) . (5) 
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FIG. 2: Time evolution of the correlation function in deposition of disks with diameters shrinking 
according to Eq. ([2]). The dashed line corresponds to t = r. The solid lines with increasing peak 
values at r/D{od) approaching 1, correspond, respectively, to t = 2r, 4r, 6r, 8t, 12t, 18t, 24r. 



At finite times the correlation function has a peak at a value r/D(oo) > 1, and actually 
vanishes at the disk contact, when r = D(t). This "correlation hole" behavior is in contrast 
with the ordinary RSA, e.g., Q], {m], for which the correlation function instead of developing 
a peak close to contact, actually increases as the disk separation decreases towards contact, 
and reaches, for finite times, a constant, non-zero value at contact. The "correlation hole" 
property represents a formation of a certain degree of short-range ordering in the system, 
and specifically, avoidance of particle clumping. 

In the limit t — > oo, the correlation hole closes and the correlation function seems to 
develop a weak divergence asr-> D(oo). In fact, for constant disk diameters, this property 
has been studied by asymptotic analytical arguments 44] and numerically [51]: The 
weak logarithmic divergence of the correlation function of RSA at contact, is not easy to 
quantify numerically, and our data were not accurate enough to study this limit for shrinking- 
diameter disk deposition by directly estimating P2{r, 00). However, the time-dependence of 
the peak values in Figure [2J follows a logarithmic divergence as t increases, as shown in 
Figure [3j This is reminiscent of a similar logarithmic divergence predicted 44J for the values 
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FIG. 3: Time-dependence of the peak values of the correlation function, demonstrating the expected 
logarithmic divergence as t — > oo. The least-squares fit of the peak values, shown here for times 
t = 4r, 6t, 8t, lOr, 12r, 18r, 24r, yields the slope of the solid line as 0.37 ± 0.02 (but note a 
comment in Subsection IIII Cl below regarding low reliability of this "error-bar" estimate). 



of P 2 (D,t) at contact (r = D) for fixed-D RSA. 

The formation of the finite-time correlation hole at disk-contact, was observed also for the 
non-exponential disk-diameter time-dependences denned in Eqs. (J3JSJ) ; see Figure HI This 
property of the correlation function is qualitatively similar for all three time-dependence 
protocols studied. However, we did not study larger time values for the non-exponential 
time-dependences. 

III. DEPOSITION OF SEGMENTS ON A LINE 
A. Numerical results 

In this section we report numerical results for the ID deposition of segments on an initially 
empty line. The segment length is a function of time, £(t), monotonically decreasing from 
£(0) to £(oo). Our numerical results in ID were obtained for the exponential time dependence 
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FIG. 4: The correlation function for t = lOr, in deposition of disks with diameters shrinking 
according to Eq. ([2]): dashed line; Eq. ©: solid line; and Eq. dotted line. 



similar to the 2D case, 

£(t) =^(oo){l + exp[-i^(oo)t]}. (6) 

We denote the number of deposition attempts per unit time per unit length (the flux) by 
R. The arriving segments are adsorbed only if they do not overlap any previously deposited 
ones; see Figure [TJ The quantity of interest is the gap density distribution function, G(x,t): 
The density of gaps (measured between the ends of the nearest-neighbor deposited segments) 
of length between x and x + dx at time t is G(x, t)dx. The density of deposited segments at 
time t, which in ID equals the density of gaps, is given by 

POO 

n(t) = / G(x,t)dx. (7) 
Jo 

Our numerical simulations followed the same procedure as in 2D. In ID, we used periodic 
boundary conditions and system sizes 10000£(oo). The time scale, cf. Eq. fl5J), was redefined 
for ID, 

t = l/M(oo) . (8) 

Our results, shown in Figure El represent averages over 10000 independent Monte Carlo 
runs. One can see that the gap density distribution function is a smooth function of x for 
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FIG. 5: Time evolution of the dimensionless gap-size distribution function £ 2 (oo)G(x, t) for seg- 
ments with length shrinking according to Eq. ©. The dashed line corresponds to t = r. The solid 
lines with increasing peak values at x/£(oo) eventually approaching 0, correspond, respectively, to 
t = 2r, 4r, 6r, 8r, 12r, 18r, 24r. 



finite times, t, and for larger times it develops a pronounced peak at x < £(t) approaching 
x = 0. However, at x = for finite times there is a "correlation hole." In the large-time 
limit, G(x, t — > oo) appears to develop a step at x = £(oo), as well as a divergence at x = 0. 
These properties will be further discussed in the next two subsections. 



B. Analytical considerations 

In this subsection, we will consider an analytical kinetic equation approach to the ID 
deposition problem. This will allow us to elucidate the origin of the "correlation hole" 
at x = 0. Kinetic equations, closed- form or hierarchies of, for RSA can be formulated 
for the gap density distribution function as well as for several other correlation-function- 



type quantities. They have been widely used 



approximation schemes; see reviews 



in Q, Q, 3, Q. 



or deriving exact results and low-density 



For the present problem, we have 



-R9(x-£(t))[x-£(t)]G(x,t) + 2R r G{y,t)dy + £(t) dG [ X,t ^ . (9) 



9t Jx+£(t) ^X 
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Here the first term represents the destruction of gaps of length x by segment deposition, with 
#(...) denoting the standard step function which vanishes for negative arguments and is 1 
for positive arguments. The second term represents the creation of gaps of length x due to 
deposition of segments landing in larger gaps of size y > x + £(t). The last term describes the 
increasing size of all the g result of segment shrinkage. Note that the case £(t) = is 

exactly solvable 53|, and we will be using some of the exact results in the next subsection. 
For large x, we obviously have the boundary condition G(x — > oo,t > 0) = for Eq. (Q. 
The following discussion explores the form of the boundary condition at the small- £ side of 
the distribution. 

There are two kinetic processes that alter gap lengths in our problem. First, deposition 
of segments can reduce the length by £(t). This process can generate very small, positive 
gap lengths, e, from available gaps of sizes £(t) + e. At the same time, the second process, 
that of segment shrinkage, causes all gap lengths to "drift" towards larger positive values 
with velocity —£{t). These two kinetic processes determine the evolution of the gap-length 
distribution G(x,t) from its initial value G(x, 0). While our numerical simulations were 
for G(x, 0) = 0, we note that in principle one can generate other translationally invariant 
initial distributions by an appropriate preparation of the initial state. In fact, one can even 
prepare initial distributions that extend to small negative gap lengths. This entails allowing 
overlaps for the segments initially placed on the line, but not for those later depositing. If 
the segments are ordered according to their center-point positions, and each gap is measured 
as a consecutive center-point distance minus £{t), then the initial-distribution overlaps, even 
if some are multi-segment, can be unambiguously counted as negative gaps. The two kinetic 
processes will then remain the same: Gaps larger than £(t) can be shortened due to deposition 
events, while at the same time all the gaps (positive, zero, negative) also increase towards 
positive values, due to segment shrinkage. 

The above considerations suggest that equation Eq. ([9]) strictly speaking should be math- 
ematically considered for — oo < x < oo. An attempt to limit it to < x < oo, and also 
use moment definitions, such as the zeroth-order moment Eq. (I7j), with integration over 
< x < oo, may in general yield wrong results: Neglecting the "flow of length" from the 
negative-x values may violate length conservation. However, in our case, specifically for the 
initial condition G(x, 0) = (and with nonvanishing £{t) for all t, including in the t — > oo 
limit, which avoids a possible singular limit), it is obvious that the problem should be de- 
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finable with a boundary condition at x = because all the kinetics of the process occurs 
in < x < oo. Indeed, together with the zeroth-order moment, the first-order moment of 
the distribution can be used to determine the applicable boundary condition directly from 
lengt h- conservation , 

POO 

£(t)n{t) + xG(x,t)dx = l. (10) 
Jo 

Here the first term is the density of the covered area, whereas the second term is the density 
of the uncovered area. Both densities are "per unit length" in ID, and therefore they sum 
up to 1. Differentiation with respect to time, t, and the use of Eq. (jUJ), then yield, after 
some algebra, the result 

£(t)G(0,t) = . (11) 
For shrinking segments, i < 0, we are thus led to our main conclusion, 

G(0,t) = 0, (12) 

except perhaps in the limit t — > oo (in which £ vanishes). This result indicates that the 
"correlation hole" for finite times, at x = is generic for the shrinking- segment RSA, no 
matter how fast is the deposition kinetics that tends to "fill the hole" in the distribution 
at small x: The presence of the correlation hole (the depletion of the gap-distribution near 
x = 0) follows from the boundary condition established for x = 0. 



C. Properties of the gap-length distribution 

For ordinary ID RSA of segments of fixed length £(t) = I = ^(oo), Figure [6] illustrates the 
exact solution {53I for several times and in the limit t — > 00. Specifically, the gap distribution 
has a finite value at x = 0, consistent with Eq. ( ITTj) for I = 0. This value actually diverges 
in the t — > 00 limit, while a logarithmic singularity develops near x = 0. This suggests that 
ordinary RSA is actually a rather singular limit of the more general shrinking-segment RSA. 
Indeed, the point x = £(t) corresponds to discontinuity in the first, deposition term in the 
kinetic equation, Eq. (J9J). As a result, the fixed-£ RSA gap-distribution has a discontinuous 
derivative at x = £, which becomes an actual discontinuity (jump) in the t —> 00 limit. On 
the other hand, with segment shrinking allowed the distribution is apparently continuous 
and smooth (as numerically observed), likely analytic at all points internal to the domain 
of definition, < x < 00, of the gap distribution. Adding the segment shrinkage process 
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FIG. 6: Time evolution of the dimensionless gap-size distribution function £ 2 G(x, t) for segments 
of constant length t. The dashed line corresponds to t = r. The solid lines correspond to t = 2r, 
4t, 6r, 8t, with increasing peak values at x = 0, increasing values at x = £, and with respectively 
smaller values for x 3> £■ The dotted line corresponds to t — > oo, with the values for x > £ equal 
in this limit. 



seems to smooth the singularities out for finite times; cf. Figure [5] In terms of the kinetic 
equation, this is a consequence of that the added third, ~ dG/dx term plays the role of 
the "diffusive" (second-derivative) smoothing contribution with respect to the ^-dependent 
integral (the second term). This diffusive property of the kinetic equation suggests that 
the function G(x,t) has no singularities for any x > for finite time, though in the limit 
t — > oo (when £{t) — > 0) the divergence as x — > + and the discontinuity at x — £(oo) are 
asymptotically restored. 



The exact solution technique for the constant-length case 



531 ] involves the use of the 



exponential-in-x ansatz for G(x, t) for x > £, recently detailed in applications for related 



models in 



271 ]. (281 ] . and then solution of the x < £ equation, Eq. Q, which becomes tractable 



because the first term is not present, whereas the integration in the second term involves 
a simple exponential integrand. Thus, this approach by its nature yields discontinuities if 
not in the function then in its first or possibly higher-order derivatives. Attempts to use 
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FIG. 7: The dashed straight line is the asymptotic large-time form of the gap-distribution function 
for the fixed-segment-length RSA; its slope is 2e~ 2 " 1 see Eq. (|13p . The solid straight line is a fit of 
our data (shown as symbols) for t > 4r, for shrinking-segment RSA, with the estimated slope of 
0.51 ±0.02. 



this approach, as well as more complicated, exponential-multiplied-by-polynomial ansaetze 
for x > £(t), for time- dependent £(t) yield solutions which satisfy the equation but possess 
unphysical discontinuities. We consider it unlikely that the shrinking-segment RSA problem 
in ID can be solved exactly by the presently known techniques. 

Finally, let us consider the extent to which our numerical data can confirm the restoration 
of the logarithmic singularity near x = in the t — ► oo limit. As in 2D, our data are 
not accurate enough to confirm the form of the x-dependence of the expected logarithmic 
divergence for x <C £(oo), i.e., G(x,oo) oc —]n(x/£(oo)). Regarding the t-dependence, 
Figure [7] shows our ID data for the peak values at small positive x, see Figure O as well 
as the asymptotic large-time form of the exact solution peak values at x = for constant 



(oo), the latter given by the relation 53 ] 



£ 2 G(t,0) ~ 2e- 27 ln(t/r), £>r, (13) 

where 7 = 0.5772156649 ... is Euler's gamma constant. Typical for data fits for logarithmic 
dependences, the estimated slope given in Figure may gradually drift as the t-variable 
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varies over many decades. It is possible that it reaches the larger value of the fixed-segment 
RSA: We consider the error estimate, 0.02, here and in the 2D case (Figure [3]) inconclusive. 

In summary, we investigated an RSA model in which deposition of arriving objects, which 
tends to form small gaps/voids, competes with the process of shrinking of the objects, both 
those already deposited and those newly arriving. The latter process always "wins" at finite 
time, and as a result "correlation hole" opens up in the correlation properties considered. 
In fact, analytical considerations in ID suggest that this property is generic for an initially 
empty substrate. Numerical Monte Carlo results reported, support this conclusion for the 
time-dependences studied. For the exponential- decay time dependence of the object sizes 
studied in ID and 2D, we also found preliminary evidence that logarithmic divergence of the 
correlation properties at contact is restored as t — > oo. For ID, the discontinuous behavior 
of the gap-distribution at x — £(oo) also reemerges in this limit. 

The authors thank Drs. S. Minko and I. Tokarev for instructive discussions, and acknowl- 
edge support of this research by the NSF under grant DMR-0509104. 
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